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ABSTRACT 

We report spectroscopic observations of the star HIP 50796, previously con- 
sidered (but later rejected) as a candidate member of the TW Hya association. 
Our measurements reveal it to be a single-lined binary with an orbital period of 
570 days and an eccentricity of e = 0.61. The astrometric signature of this orbit 
was previously detected by the HIPPARCOS satellite in the form of curvature in 
the proper motion components, although the period was unknown at the time. 
By combining our radial velocity measurements with the HIPPARCOS interme- 
diate data (abscissae residuals) we are able to derive the full three-dimensional 
orbit, and determine the dynamical mass of the unseen companion as well as a 
revised trigonometric parallax that accounts for the orbital motion. Given our 
primary mass estimate of 0.73 Mq (mid-K dwarf), the companion mass is deter- 
mined to be 0.89 Mq, or ~20% larger than the primary. The likely explanation 
for the larger mass without any apparent contribution to the light is that the 
companion is itself a closer binary composed of M dwarfs. The near-infrared 
excess and X-ray emission displayed by HIP 50796 support this. Our photomet- 
ric modeling of the excess leads to a lower limit to the mass ratio of the close 
binary of g ~ 0.8, and individual masses of 0.44 — 0.48 M© and 0.41 — 0.44 M©. 
The new parallax (tt = 20.6 ±1.9 mas) is significantly smaller than the original 
HIPPARCOS value, and more precise. 

Subject headings: binaries: spectroscopic — binaries: visual — stars: individual 
(HIP 50796) — stars: late-type — techniques: radial velocities — techniques: 
spectroscopic 
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1. Introduction 

In recent years there has been considerable interest in the study of nearby groupings of 
young stars that do not appear to be associated with any molecular clouds. The prototypical 
example is perhaps the TW Hydrae association (Kastner et al. 1997), with an estimated age 
of 8-10 Myr and some two dozen members. Searches for new members of this group have 
relied on X-ray properties, kinematics, infrared (2MASS) colors, and spectral features (e.g., 
Makarov & Fabricius 2001; Gizis 2002; Tachihara et al. 2003; Song, Zuckerman & Bessell 
2003), and have produced many candidates and motivated detailed follow-up studies to 
confirm them. 

One of those searches, by Makarov & Fabricius (2001), used proper motions from the 
Tycho-2 catalog (H0g et al. 2000) for X-ray sources hsted in the ROSAT Bright Source Cat- 
alog (Voges et al. 1999). The authors considered a large area surrounding the previously 
known members of the TW Hya association, and proposed a list of 23 candidate mem- 
bers that were subsequently followed up spectroscopically by Torres, Neuhauser & Latham 
(2001), Song, Bessell & Zuckerman (2002), and Torres et al. (2003). The present paper 
discusses one of those objects, HIP 50796 (also known as BD-09 3055, GSC 05493 00324, 
IRXS J102219.2-103302, TYC 5493-324-1, a = 10'^22'"17!99, 5 = -10°32'15'.'5, J2000, 
V — 10.80, SpT = K). The kinematic model of the association constructed by Makarov & 
Fabricius (2001) made specific predictions as to the radial velocity expected for each candi- 
date. Although the measurements for HIP 50796 by Torres, Neuhauser & Latham (2001) 
happened to agree perfectly with those predictions (+13.1 ± 1.0 kms~^). Song, Bessell & 
Zuckerman (2002) found that the star shows no detectable Li A6708 absorption in its spec- 
trum (equivalent width < 10 mA) , which indicates it is not very young and therefore essen- 
tially rules it out as a true member. They also measured a radial velocity of +22.4±0.9 km s~^ 
that is nearly 10 kms~^ different from the previous measurements, and proposed that the 
object is a spectroscopic binary. Further measurements by Torres et al. (2003) confirmed 
the velocity variations and indicated an orbital period of perhaps 1 or 2 years. 

Even though it was clear that HIP 50796 is not related the TW Hya association, the 
object was kept on the observing hst at the Harvard-Smithsonian Center for Astrophysics 
(CfA) for the purpose of establishing the orbit. Preliminary solutions appeared to imply a 
mass for the companion that is larger than the primary, which was puzzling given that the 
spectrum is single-lined. Additionally, the HIPPARCOS satellite (ESA 1997) detected cur- 
vature in the proper motion further supporting the binary nature of the object. We present 
here the complete analysis of HIP 50796 that combines the astrometry and the spectroscopy 
and allows us to investigate in more detail the nature of the overmassive secondary star. 
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2. Spectroscopic observations and reductions 

HIP 50796 was observed spectroscopically at the CfA with an echelle spectrograph on 
the 1.5-m Tilhnghast reflector at the F. L. Whipple observatory on Mt. Hopkins (Arizona). 
A single echelle order was recorded using a photon-counting intensified Reticon detector at 
a central wavelength of 5187 A, giving a spectral coverage of 45 A. The resolving power is 
A/AA 35,000. A total of 34 spectra were obtained from April 2002 until May 2005, and 
the signal-to- noise (S/N) ratios achieved range from 5 to about 30 per resolution element 
of 8.5 kms"-*". Four archival spectra (S/N of 9-30) were also used, and were obtained for 
another program in 1985-1986 with nearly identical instrumentation on the Multiple Mirror 
Telescope (also on Mt. Hopkins, Arizona), prior to its conversion to a monolithic mirror. 

Radial velocities were derived using XCSAO (Kurtz & Mink 1998), a one-dimensional 

cross-correlation program well suited to our relatively low S/N spectra that runs under 
IRAF^. For the template we used a synthetic spectrum selected from an extensive library of 
spectra based on model atmospheres by R. L. Kurucz^ (see Nordstrom et al. 1994; Latham 
et al. 2002). These calculated spectra are available for a wide range of effective temperatures 
(Teff), projected rotational velocities (f sini), surface gravities (logs') ^^'^ metallicities. The 
optimum template was determined from grids of cross-correlations over broad ranges in 
Teflf and vsini, seeking to maximize the average correlation weighted by the strength of each 
exposure. Solar metallicity was assumed to begin with. The above calculations were repeated 
for a range of logs' values from 2.0 to 5.0, which allowed us to derive a rough estimate of the 
surface gravity for HIP 50796. We obtained logs = 4.2 ± 0.5 and %s = 4600 ± 150 K. The 
entire procedure was repeated for metallicities different from solar, but the best match was 
achieved for templates with the solar composition. The temperature we derive corresponds 
to a main-sequence star with spectral type mid-K (Cox 2000). The rotational broadening 
was found to be neghgible (formally vsini — 1 ± 3 kms~^). Song, Zuckerman & Bessell 
(2003) reported a considerably larger value of 8 kms~^. 

The stability of the zero-point of our velocity system was monitored by means of expo- 
sures of the dusk and dawn sky, and small run-to-run corrections were applied in the manner 
described by Latham (1992). The radial velocities in the hehocentric frame including these 
corrections are listed in Table 1. The typical precision of a single measurement is 0.5 kms~^. 
Variations over the 19.4 years of coverage are obvious and show a peak-to-peak amplitude of 
about 40 kms~^. The period is 570 days, and the eccentricity (e = 0.61) is quite significant. 



-'^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associ- 
ation of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. 

^Available at http://cfaku5.cfa.harvard.edu. 
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An orbital solution based on these velocities is presented in the second column of Table 2. 
Although it is no longer relevant in connection with the possible membership to the TW Hya 
association, we mention in passing that the ccnter-of-mass velocity 7 = +24.77 ±0.14 kms~^ 
of HIP 50796 is far from the value of +13.1 kms^^ that had been predicted by Makarov & 
Fabricius (2001). The most significant finding, however, is the large mass function of the 
binary: f{M) = 0.258 + 0.011. For a typical primary mass corresponding to the temperature 
we derive, the mass function implies a secondary that is more massive than the primary, yet 
no signs of another star are obvious in our spectra. We discuss this further below. 

3. Astrometric observations 

HIP 50796 was observed by the HIPPARCOS satellite from December 1989 to November 
1992, during which 53 astrometric measurements were made. The trigonometric parallax 
was determined to be tthip = 29.40 ± 2.69 mas. The uncertainty is somewhat larger than 
usual, but perhaps understandably so given that the star is relatively faint {V — 10.80). 
Additionally, the HIPPARCOS reductions revealed significant curvature in the proper motion 
in the amount oid^a/dt = —10.90+6.51 mas yr~^ and djis/dt — +22.33+5.84 mas yr~^, most 
likely a reflection of the orbital motion of the object. Given its period of 570 days, however, 
nearly two cycles of the orbit were covered during the HIPPARCOS campaign (which lasted 
1058 days) so that the meaning of those coefficients is obscured. The proper motions and 
the parallax as reported in the HIPPARCOS catalog are also likely to be compromised. 

Since the satellite observations are sensitive enough to have detected the orbital motion 
on the plane of the sky, the best course of action is therefore to make use of the individual 
astrometric measurements ("abscissae residuals"; ESA 1997) together with the velocities to 
constrain the orbit, and at the same time re-derive the position, proper motion, and parallax. 
Those HIPPARCOS measurements are hsted in Table 3. We describe the global solution 
below. 



4. Orbital solution 

The radial velocities allow for the determination of the usual spectroscopic elements, 
which are the period (P), center-of-mass velocity (7), semi-amplitude of the velocity variation 
{K), eccentricity (e), longitude of periastron for the primary (cui), and time of periastron 
passage (T). On the other hand, the visual elements constrained by the astrometry are the 
angular semimajor axis (a), the inclination angle (i), the position angle of the ascending 
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node (Q), as well as P, e, uji, and T. The combination of the two kinds of measurements 
thus provides redundancy in the four latter orbital elements that strengthens the solution. 
Furthermore, the fact that we see no indication of the secondary in the spectrum of HIP 50796 
means that the motion seen by HIPPARCOS should essentially be that of the primary star 
around the center of mass of the binary. This introduces a relation between the angular 
semimajor axis of the primary, ai, and its velocity amplitude, given by 

ai = 9.191967 x 10~^ • ttKPVI -e^/sini , (1) 

where ai and the trigonometric parallax vr arc expressed in milli-arc seconds, K in kms~^, 
and P in days. We use this relation to eliminate the angular semimajor axis (oi) as an 
unknown. Since the HIPPARCOS observations are made in an absolute frame of reference, 
the abscissae residuals contain information on the parallax as well as the position and proper 
motion of the barycenter of the binary. Five additional variables are thus introduced (Aa*, 
AS, A/x* , AfjLs, Att)^, which represent corrections to a fiducial point of reference that yield 
improved estimates of those quantities over the catalog values (see ESA 1997). 

A total of 13 variables were adjusted in the solution, using standard non-linear least- 
squares techniques (Press et al. 1992). The formalism used for incorporating the abscissae 
residuals from HIPPARCOS into the fit follows closely that described by van Leeuwen & 
Evans (1998) and Pourbaix & Jorissen (2000), and is described in more detail in the Ap- 
pendix. The relative weights of the spectroscopic and astromctric observations were adjusted 
to yield separate reduced values near unity. In this way we determined a scale factor for 
the original uncertainties of the HIPPARCOS abscissae residuals of 1.14, and a factor of 1.43 
for the original velocity errors. 

Due to the faintness of the target the median error of a single HIPPARCOS measurement 
is 5.7 mas, so that the constraint provided by the astrometry is relatively weak. Nevertheless, 
the results of this initial fit yielded improvements in all of the orbital elements in common 
with the spectroscopy. The value derived for the angular semimajor axis of the primary was 
Oi = 16.9±1.6 mas. As a test to see if the astrometry is able to detect the light contribution 
from the companion, we ran another solution in which we left the semimajor axis of the 
apparent orbit as a free parameter. In this case that semimajor axis (ophot) corresponds 
to the motion of the center of light of the binary, or photocenter, rather than that of the 
primary, and it should in principle be smaller if the light from the secondary is significant. 
This fit gave a value of Ophot = 18-1 ±3.6 mas, which is not only not smaller but is also much 
more uncertain (because it is no longer constrained by the spectroscopy through eq.[l]). 



^Following the practice of the HIPPARCOS catalog we define Aa* = AacosJ and A/x* = Afia cos 5. 
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Nevertheless, it is still consistent with the previous result, within the errors. We conclude 
from this that the companion is not bright enough to produce a measurable astrometric 
effect, and we proceed under the assumption that that star is invisible (but see §6). 

Slight differences were found in the proper motion components (//* , /is) compared to the 
values reported in the HIPPARCOS catalog, which was not unexpected given that the catalog 
values do not account for orbital motion. An external check on these motions is available 
from the Tycho-2 catalog (H0g et al. 2000), in which /i* and fis are derived by combining the 
HIPPARCOS position (mean epoch 1991.25) with positions from ground-based meridian- 
circle and photographic catalogs going back several decades, and up to a century in some 
cases. Given that the period of HIP 50796 is only 570 days, the orbital motion over a baseline 
of decades should tend to average out in the Tycho-2 analysis, resulting in a more accurate 
measure of the proper motion than that reported by HIPPARCOS. Indeed, as shown in 
Table 4 the values of //* and ns from our initial fit are much closer to the Tycho-2 values 
than the HIPPARCOS values. This also suggests that the solution might benefit if we made 
use of the constraint provided by Tycho-2. Consequently, we incorporated the Tycho-2 
proper motions into the fit as measurements, along with their uncertainties. The resulting 
proper motions from this combined solution are listed for comparison in the fourth row of 
Table 4. The remaining orbital elements are hardly affected, except for the slightly smaller 
uncertainties. 

The complete results of this combined fit arc given in the third column of Table 2. 
The scmimajor axis of the relative orbit, a, is inferred to be 30.8 mas (1.58 AU in linear 
units). The orbital period is determined to better than 0.1% by virtue of the 19.4-yr baseline 
provided by the velocity measurements (more than 12 cycles). Perhaps one of the most 
significant improvements coming from the combination of the astrometry and spectroscopy 
is in the trigonometric parallax. Compared to the value from the HIPPARCOS catalog, 
our parallax that takes full account of the orbital motion is nearly 10 mas smaller, which 
corresponds to a 50% change in the distance to HIP 50796 (see Table 4), and is considerably 
more precise. This has important consequences for the luminosity estimates discussed below. 

According to our solution the orbit of the binary is seen nearly edge-on: the inclination 
angle is 85° ± 13°. Although in principle this would allow for eclipses, these are very unlikely 
given the 570-day period. Inspection of the HIPPARCOS photometry for HIP 50796 shows 
no evidence of any eclipse events (four conjunctions occurred during the mission), and the 
peak-to-peak brightness variation recorded is ~0.1 mag, which is probably at the level of the 
uncertainties for a star as faint a.sV = 10.80. Nevertheless, carefully scheduled observations 
may perhaps be of interest since the satellite measurements do not have the optimal sampling. 
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5. The nature of the companion 

The information provided by our combined solution allows a direct determination of 
the mass of the unseen companion of HIP 50796, given an estimate of the primary mass. 
For this we use our effective temperature determination in §2 (Tgff = 4600 ± 150 K), and 
the brightness of the object. The latter was measured during the HIPPARCOS mission as 
Hp — 10.9366 ± 0.0055, in the passband of the satellite. Conversion to the Johnson system 
(ESA 1997) yields V — 10.80 ±0.02, where the uncertainty is our conservative estimate that 
includes the transformation. With our parallax of tt = 19.5 ± 1.8 mas (corresponding to a 
distance of 51.3 ± 4.9 pc) the absolute magnitude in the visual band is My = 7.25 ± 0.20, 
ignoring extinction, in which essentially all of the uncertainty comes from the parallax error. 

The location of the primary star in the H-R diagram is shown in Figure 1. Model 
isochrones from the series by Girardi et al. (2000) for solar metallicity agree well with the 
measured properties of HIP 50796, quite independently of age (ranging from 1 Gyr to 5 Gyr 
in the figure). From these models we estimate the primary mass to be Mi = 0.73 ±0.05 Mq, 
in which the uncertainty accounts for observational errors as well as age. The surface gravity 
inferred from the models is logg = 4.65 ± 0.05, in agreement with our crude estimate 
from §2. Although the precise metallicity of the star is unknown, our tests in §2 suggested 
a composition near solar. Additionally, the space motion of HIP 50796 in the Galactic 
frame based on our parallax, proper motion, and center-of-mass velocity is quite small {U — 
—8 kms~"^, V = —14 kms~^, W = ±9 kms~^, relative to the Local Standard of Rest), which 
supports a Population I origin and a heavy element abundance probably close to that of the 
Sun, as we have assumed. 

The adopted mass for the primary leads to a companion mass of Mpomp = 0.88±0.05 Mq. 
The uncertainty in the primary mass contributes only about 25% to this error. The secondary 
star is thus 20% more massive than the primary, which implies it cannot be a single main- 
sequence star. A more careful examination of our spectra was carried out using the two- 
dimensional cross-correlation algorithm TODCOR (Zucker & Mazeh 1994), to place limits on 
the brightness of the companion. This technique often allows the detection and measurement 
of faint secondaries that are difficult to see in standard one-dimensional cross-correlation 
diagrams. No evidence of another set of lines was found down to the level of the noise, at 
about 10% of the brightness of the primary. Additionally, Hemenway et al. (1997) reported an 
observation of HIP 50796 by speckle interferometry in which no companions were apparently 
seen in the range from 30 mas to 0'/7. Details of this measurement are unavailable. 

One possibility is that the companion is a white dwarf. Depending on its temperature, 
it may be possible to detect it in the ultraviolet, although no such observations are available 
to our knowledge. Given the present mass of the object, initial-final mass relations for white 
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dwarfs (e.g., Ferrario et al. 2005, and references therein) suggest a progenitor of approxi- 
mately 4-5 Mq, corresponding to a mid-B star. Alternatively, the companion of HIP 50796 
could itself be a (closer) binary composed of lower main-sequence stars, making the system 
a hierarchical triple. In this case, however, the mass ratio of the binary cannot be too small 
or we would have seen the brighter of the two stars in our spectra''. 

Possible signs of this type of configuration might be seen in the form of an infrared 
excess. JHKg measurements for HIP 50796 are available from the 2MASS catalog. With 
the visual magnitude listed earlier, and after conversion to the standard Johnson system as 
defined by Bessell & Brett (1988) using transformations by Carpenter (2001), wc obtain the 
observed colors V — J, V — H, and V— K listed in the first line of Table 5. The predicted 
colors for a single star of the assumed mass (Mi = 0.73 M©) according to the models by 
Girardi et al. (2000) for a representative age of 2 Gyr are indeed bluer, by as much as 
0.4 mag. However, there are a number of reasons to doubt the theoretical colors in this case, 
including indications of missing opacity sources and other limitations in the physics (e.g., 
Baraffe et al. 1998; Delfosse et al. 2000; Chabrier et al. 2005, and references therein), and 
especially disagreements with empirical mass-luminosity relations, particularly in the visual 
band. Therefore, we have chosen here to rely on the empirical relations. Those by Henry 
& McCarthy (1993) give redder colors than Girardi et al. (2000) for the primary (see third 
line of Table 5), but are still considerably bluer than observed, by 0.17 mag, 0.29 mag, and 
0.25 mag loxV—J,V— H, and V— K, respectively. These differences are 6-7 times larger 
than the observational errors, suggesting the infrared excess may be real. 

With the aid of the Henry & McCarthy (1993) relations, we have modeled this excess 
by computing the visual and JHK magnitudes of the components of the unseen companion 
(hereafter star 2 and star 3) over a range of possible masses for each star, subject to the 
constraint that their sum be the value we determined above, Mpomp = 0.88 Mq. We then 
added the light of these stars to that of the primary (star 1) in each passband, and calculated 
the combined colors for the triple system. For convenience we parameterized the problem 
in terms of the mass ratio of the binary (g = M3/M2, where star 2 is the more massive 
component). Additionally, we allowed for changes in the mass of the main star, and sought 
to produce the best match to the observed colors as well as the combined absolute visual 
magnitude, in a sense. A solution was found for a mass Mi = 0.728 Mq and a binary 
mass ratio q — 0.795, which reproduces the observed visual brightness and all infrared colors 
within the uncertainties (see last row of Table 5). Contours of the surface as a function of 



A possible exception would be if the orbital period of the binary is short enough that the stars are spun 

up by tidal forces and arc rotating very rapidly, in synchronism with the orbital motion. In that case the 
spectral lines may be broad enough to reduce the contrast and escape detection. 
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Mi and q are shown in Figure 2, where the best fit is indicated by a dot. The masses of the 
binary components from this model arc M2 = 0.49 Mq and M3 = 0.39 Mq. The predicted 
magnitude difference in the visual band between the brighter of these stars and star 1 is 
3.0 mag, which is below our threshold for spectroscopic detection (~2.5 mag; see above) and 
is thus consistent with all observational constraints. Although this is formally the best fit 
to the observations, we note that the surface is nearly fiat in the q direction for q larger 
than about 0.7-0.8 at a fixed value of Mi, so that larger mass ratios up to and including 
unity are equally acceptable and reproduce the observed colors just as well (see Figure 2). 
The mass of star 1, on the other hand, is very tightly constrained. At g = 1 stars 2 and 3 
would be at their faintest, each about 3.3 mag fainter than star 1. 



6. Constrained orbital solution 

Even though, according to the model above, both stars in the binary are too faint to 
be seen individually in our spectra, their combined light is not negligible compared to the 
brightness of star 1. The difference is approximately 2.5 mag in the visual band. This 
suggests that our earher assumption that the astrometric motion detected by HIPPARCOS 
traces only the motion of star 1 may not be entirely correct. The extra light (even with 
AV = 2.5) must affect the astrometry at some level, despite our failure to detect the effect 
directly, as described in §4. This will cause the semimajor axis of the photocenter, flphot, 
to be slightly smaller than that of the primary alone by a factor [B — j3)/ B ^ 0.17, where 
B — Mcomp/(Mi + Mcomp) IS the mass fraction of the companion and (3 — [1 + io° '^^^)~^ 
is its fractional fight. We therefore adjusted our orbital solution to allow for this effect, 
and recomputed all orbital elements. The results of this "constrained" solution are seen 
in the last column of Table 2. The semimajor axis of the photocenter is reduced from 
Ophot — di — 16.9 ±1.5 mas to Cphot = 14.9 ±1.3 mas, a small effect. 

The residuals of the radial- velocity observations from this fit are fisted in Table 1, and 
those of the HIPPARCOS measurements are given in Table 3. Figure 3 displays the spec- 
troscopic measurements along with the velocity curve as a function of orbital phase. The 
astrometric observations on the plane of the sky are illustrated in Figure 4 and Figure 5, in 
which the axes are parallel to the Right Ascension and Declination directions. The curious 
pattern in Figure 4 is the result of the combined effects of annual parallax, proper motion, and 
orbital motion. The dominant contribution is from the proper motion (80 mas yr^^), which 
is indicated with an arrow. Parallax and orbital motion are smaller and comparable effects. 
In Figure 5 we have subtracted the proper motion and parallactic contributions, leaving only 
the orbital motion with the 570-day period and a semimajor axis of Cphot — 14.9 mas. The 
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direction of motion is direct (counterclockwise). The individual HIPPARCOS observations 
are represented schematically in both of these figures, but are seen more clearly in Figure 5. 
Because they are one-dimensional in nature (ESA 1997), the exact location of each measure- 
ment on the plane of the sky cannot be shown graphically. The filled circles represent the 
predicted location on the computed orbit. The dotted lines connecting to each filled circle 
indicate the scanning direction of the HIPPARCOS satellite for each measurement, and show 
which side of the orbit the residual is on. The short hue segments at the end of and per- 
pendicular to the dotted lines indicate the direction along which the actual observation lies, 
although the precise location is undetermined. Occasionally more than one measurement 
was taken along the same scanning direction, in which case two or more short line segments 
appear on the same dotted lines. 

Compared to our previous solution, the parallax is increased slightly to tt = 20.6 ± 
1.9 mas, which corresponds to a distance of 48.5 ± 4.5 pc. The proper motion components 
change very little (Table 4). The binary mass derived with Mi = 0.73 Mq is marginally 
larger than before (Mcomp = 0.89 ± 0.06 Mq). The change in the parallax leads to a revised 
absolute magnitude for the system of My = 7.37 ± 0.20. With these updated values of 
^comp and My we repeated the photometric modeling using the Henry & McCarthy (1993) 
mass-luminosity relations, and obtained a solution not very different from the previous one 
(the same value of Mi, and q = 0.842, implying M2 = 0.48 M© and M3 = 0.41 Mq). The 
brightness difference between star 1 and the sum of stars 2 and 3 is unchanged, making 
another iteration of our orbital solution unnecessary. 



7. Discussion and concluding remarks 

HIP 50796 is an interesting illustration of the complementarity of spectroscopic and 
astrometric observations, and in particular of the utihty of the HIPPARCOS intermediate 
data (see also Pourbaix 2004; Jorissen, Jancart & Pourbaix 2004; Pourbaix, Jancart & Joris- 
sen 2004). While the radial- velocity measurements clearly reveal this object to be a binary 
and provide the (single-lined) spectroscopic orbit with a period of 570 days, the combina- 
tion with the astrometry has allowed us to derive the dynamical mass of the companion 
(Afconip = 0.89 M0) and also to correct the seriously biased parallax value from the original 
HIPPARCOS reductions. No indications of the companion are detected in our spectra, yet 
the orbital solution shows that it is clearly more massive than the visible star, which is a 
normal K dwarf with M a; 0.73 Mq. While we cannot completely rule out that the secondary 
is a massive white dwarf, all the evidence points quite convincingly toward the conclusion 
that the companion is itself a closer binary composed of M dwarfs, making the system a 
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hierarchical triple. 

The infrared excess that might be expected from such a configuration appears indeed to 
be present, as indicated by our modehng of the measured visual and near-infrared (2MASS) 
photometry using empirical mass- luminosity relations by Henry & McCarthy (1993), along 
with the constraint on the total mass of the unseen companion. This modeling cannot 
determine the precise mass ratio of the close binary, but is able to place a lower limit of 
about q ~ 0.8. Smaller values would make one of the stars bright enough that it would be 
seen in our spectra. We infer masses for these stars of approximately M2 = 0.44 — 0.48 Mq 
and M3 = 0.41 — 0.44 M0. The period of the close binary is unknown. We note also that 
HIP 50796 is listed as an X-ray source in the ROSAT catalog, with \ogLx = 28.62 (Micela, 
Favata & Sciortino 1997) and log(Lx/i^boi) = -2.61 (Makarov & Fabricius 2001). This X- 
ray emission might arise naturally from the presence of the M dwarfs, which are frequently 
active, particularly if they are in a short-period configuration so that tidal forces compel the 
stars to be in synchronous rotation with the orbital motion. 

Finally, the relative orbit of HIP 50796 and its unseen companion has an angular semi- 
major axis of about 33 mas. Given the eccentricity of the orbit, separations up to about 
50 mas are possible at times, which may permit a direct detection of the secondary with 
high-resolution techniques. Alternatively, infrared spectroscopy should be able to reveal the 
presence of at least the brighter of the M dwarfs directly, given the more favorable contrast 
with the main star at those wavelengths. 

The author is grateful to P. Berlind, M. Calkins, D. W. Latham, and R. P. Stefanik for 
obtaining the spectroscopic observations used in this work, and to R. J. Davis for maintain- 
ing the CfA echelle database. B. Mason and P. Hemenway are thanked for information on 
speckle measurements. This paper benefited also from helpful comments by an anonymous 
referee. Partial support for this work from NSF grant AST-0406183 and NASA's MASSIF 
SIM Key Project (BLF57-04) is acknowledged. This research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, France, of NASA's Astrophysics Data System 
Abstract Service, and of data products from the Two Micron All Sky Survey, which is a 
joint project of the University of Massachusetts and the Infrared Processing and Analysis 
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A. Incorporating HIPPARCOS data into the global orbital solution 

The intermediate data provided with the HIPPARCOS catalog are the "abscissae residu- 
als". At;, which are the difference between the satellite measurements (abscissae) along great 
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circles and the abscissae computed from the 5 standard astrometric parameters. The stan- 
dard parameters are the position of the object (oq, Sq) at the reference epoch to = 1991.25, 
the proper motion components (/i*, fis), and the parallax (tt). We follow here the notation 
in the HIPPARCOS catalog and define = aQCOs6 and yU* = fi^cosd, to incorporate the 
projection factors. The goal of an orbital solution making use of the HIPPARCOS data is to 
reduce the abscissae residuals below the values obtained from the 5-parameter solution by 
taking the orbital motion into account^. Following Pourbaix & Jorissen (2000) the sum 
for the abscissae residuals is = H* V"^ S, where 



k-- 



^ dp. 



and H* is the transpose of H. In this expression A.v is the array of N abscissae residuals 
provided by HIPPARCOS, and dv/dpk is the array of partial derivatives of the abscissae 
with respect to the fc-th fitted parameter. The number M of parameters fitted to the 
astrometry in the general case is 12: the 5 standard HIPPARCOS parameters {pi = a^, 
P2 = So, P3 = f^a, P4 = /^<5) P5 = ^) ^-ud 7 Orbital elements (ai, P, e, i, Ui, Q, T, represented 
as pk with k = 6, ...,12). V~-^ is the inverse of the covariance matrix of the observations, 
containing the abscissae uncertainties and correlation coefficients (ESA 1997, Vol. 3, eqs. 
17-10 and 17-11) also provided with the HIPPARCOS catalog. Correlations arise because 
the same original data were reduced independently by two data reduction consortia (NDAC 
and FAST; see ESA 1997), and both results are included in the solution. 

The partial derivatives dv/dpk for A; = 1 to 5 are given in the HIPPARCOS catalog 
along with the abscissae residuals. The remaining derivatives can be expressed in terms of 
the partial derivatives of v with respect to ctg and 5q. These are (ESA 1997, Vol. 3, eq. 
17-15) 

dv dv dv dri , „ / » 

l^ = 7^*l^ + 7^7r- ' ^ = 6, ...,12, A2 

dpk da*Q dpk ddo dpk 

in which ^ and r] are in our case the rectangular coordinates of the photocenter relative to 
the center of mass of the binary on the plane tangent to the sky at {a^, 5q), which are given 

by 

i = aS + /i:(t-to)+7rP« + AX (A3) 
V = 5o + Ai5(^-^o)+vrP5 + Ay . (A4) 



^The 5 standard parameters for HIP 50796 were actually computed as part of a 7-parameter solution 
incorporating the time derivatives of the proper motion components (since HIPPARCOS detected curvature 
in the proper motion). Nevertheless, the abscissae residuals available are always derived from the 5 standard 
parameters as listed in the catalog. 
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Pa and Ps are the parallactic factors, and the terms AX — Bx + Gy and Ay = Ax + Fy 
represent the orbital motion components, where A, B, F, and G are the classical Thiele- 
Innes constants. These depend only on the orbital elements ai, i, Ui, and fl (see, e.g., van de 
Kamp 1967). x and y are the rectangular coordinates in the unit orbit given hy x — cosE — e 
and y = yl — e^sinE, with E being the eccentric anomaly. 

As described by Pourbaix & Jorissen (2000), the nature of the orbital solution is such 
that only the derivative corresponding to the semimajor axis in eq.(A2) needs to be consid- 
ered. The expression for H in eq.(Al) then reduces to 



That the semimajor axis ai has actually been eliminated as a formal adjustable parameter 
in our case is irrelevant, since it still appears in the Thiele-Innes constant but is computed 
from other elements using eq.(l). 

Finally, the for the global solution that combines astrometry and spectroscopy is 
computed by adding the term corresponding to the radial velocities, in the usual manner. 




(A5) 
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8000 7000 6000 5000 4000 

Effective Temperature (K) 

Fig. 1. — Location of HIP 50796 in the H-R diagram, against the backdrop of model 
isochrones by Girardi et al. (2000) for solar metaUicity and ages ranging from 1 Gyr to 
5Gyr. 
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Fig. 2. — Morphology of the surface resulting from our modeling of the observed V and 
JHK photometry of HIP 50796 with the sum of three stars, (a) Contours of equal x^- 
The dot represents the best fit with Mi = 0.728 M© and a binary mass ratio of g = 0.795, 
corresponding to M2 — 0.49 M0 and M3 = 0.39 Mq (see text); (b) Cross-section along the 
Ml axis, showing a well defined minimum; (c) Cross-section along the q axis, showing a very 
fiat surface in the vicinity of the minimum. Solutions at higher mass ratios and the same 
Ml are not significantly worse. 
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Fig. 3. — Radial velocity observations for HIP 50796 and fitted orbit. The dotted line 
represents the center-of-mass velocity of the binary. Errors are smaller than the size of the 
points. 
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Fig. 4. — Path of the center of hght of HIP 50796 on the plane of the sky, along with the 
HIPPARCOS observations (abscissae residuals). See text or Figure 5 for an explanation 
of the graphical representation of these one- dimensional measurements. The figure shows 
the total motion resulting from the combined effects of parallax, proper motion, and orbital 
motion according to the global solution described in the text. The arrow indicates the 
direction and magnitude of the annual proper motion. 
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Fig. 5. — Residual orbital motion of the center of light of HIP 50796 on the plane of the sky, 
after removal of the parallactic and proper-motion components. HIPPARCOS observations 
are represented as in Figure 4 (a few measurements with large residuals have been omitted 
here for clarity). Because they are one-dimensional in nature (ESA 1997), the exact location 
of each measurement on the plane of the sky cannot be shown graphically. Filled circles 
represent the predicted location on the computed orbit. The dotted lines connecting to each 
filled circle indicate the scanning direction of the HIPPARCOS satellite for each measure- 
ment, and show which side of the orbit the residual is on. The short line segments at the 
end of and perpendicular to the dotted lines indicate the direction along which the actual 
observation lies, although the precise location is undetermined. Occasionally more than one 
measurement was taken along the same scanning direction, in which case two or more short 
line segments appear on the same dotted lines. The plus sign in the figure indicates the 
center of mass of the binary, and periastron is shown with an open circle. The direction of 
motion in the orbit is direct (counterclockwise). 
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Table 1. Radial velocity measurements for HIP 50796 in the heliocentric frame. 



HJD 
(2,400,000+) 


Year 


RV 
(kms-l) 


(kms-1) 


(O-C) 
(kms-l) 


Orbital 
Phase 


46422.9750 


1985 


.9767 


+16, 


.31 


0, 


.31 


-0 


.06 





.605 


46537.6905 


1986, 


.2907 


+14, 


,05 


0, 


,57 


+1, 


.11 


0, 


,806 


46540.7560 


1986, 


.2991 


+12, 


,68 


0, 


,41 


-0, 


.25 


0, 


.812 


46569.7043 


1986, 


.3784 


+12, 


,01 


0, 


,73 


-1, 


.43 


0, 


.862 


52391.6615 


2002, 


.3180 


+47, 


.99 


0, 


.86 


-1, 


.38 





.064 


52395.7197 


2002, 


.3291 


+48, 


.10 


0, 


.59 


+0, 


.03 


0, 


.071 


52419.6668 


2002, 


.3947 


+40. 


.40 


0, 


.70 


-0, 


.96 


0, 


.113 


52421.6565 


2002 


.4002 


+40, 


.59 


0, 


.50 


-0 


.30 





.116 


52423.6634 


2002, 


.4056 


+40, 


,47 


0, 


,60 


+0, 


,05 


0, 


.120 


52424.6403 


2002, 


.4083 


+41, 


,03 


0, 


,70 


+0, 


,83 


0, 


.122 


52655.0264 


2003 


.0391 


+ 18, 


,00 


0, 


,71 


-0 


..54 





.525 


52656.9803 


2003, 


.0444 


+17, 


.85 


0, 


.50 


-0, 


.60 


0, 


.529 


52663.9288 


2003, 


.0635 


+18, 


.03 


0, 


.56 


-0, 


.09 


0, 


.541 


52691.9163 


2003. 


.1401 


+16. 


.96 


0. 


.54 


+0, 


.09 


0, 


.590 


52712.8028 


2003 


.1973 


+14, 


.63 


0, 


.97 


-1, 


.37 





.627 


52752.7176 


2003, 


.3066 


+14, 


,95 


0, 


,60 


+0, 


.44 


0, 


.697 


52772.7344 


2003, 


.3614 


+14, 


,13 


0, 


,71 


+0, 


.26 


0, 


.732 


52985.0078 


2003, 


.9425 


+43, 


.29 


0, 


,51 


+0, 


.53 


0, 


.104 


53013.9323 


2004, 


.0217 


+36, 


.91 


0, 


.44 


+0, 


.36 


0, 


.154 


53035.8873 


2004 


.0818 


+32, 


.92 


0, 


.51 


-0 


.19 





.193 


53044.9091 


2004 


.1065 


+31, 


,77 


0, 


,47 


-0 


.14 





.209 


53046.9420 


2004, 


.1121 


+31, 


,34 


0, 


,47 


-0, 


.32 


0, 


.212 


53073.7601 


2004, 


.1855 


+28, 


,56 


0, 


,44 


-0, 


.12 


0, 


.259 


53125.7468 


2004, 


.3278 


+25, 


.00 


0, 


.44 


+0, 


.67 


0, 


.350 


53132.7248 


2004, 


.3470 


+24, 


.29 


0, 


.44 


+0, 


.45 


0, 


.362 


53337.0326 


2004, 


.9063 


+13, 


.83 


0, 


.46 


-0, 


.25 


0, 


.720 


53362.0644 


2004 


.9749 


+13, 


.70 


0, 


.46 


+0, 


.32 


0, 


.764 


53371.9950 


2005, 


.0020 


+12, 


,45 


0, 


,39 


-0, 


.72 


0, 


.782 


53393.8773 


2005, 


.0620 


+13, 


,90 


0, 


,39 


+0, 


.97 


0, 


.820 


53420.9096 


2005, 


.1360 


+13, 


,88 


1, 


.13 


+0, 


,47 


0, 


.867 


53426.8509 


2005, 


.1522 


+13, 


.46 


0, 


.51 


-0, 


.27 


0, 


.878 


53428.8623 


2005, 


.1577 


+13, 


.35 


0, 


.51 


-0, 


.52 





.881 


53451.7621 


2005, 


.2204 


+17, 


.17 


0, 


.36 


+0, 


.11 


0, 


.921 


53457.7802 


2005 


.2369 


+19, 


.09 


0, 


.50 


+0, 


.37 


0, 


.932 


53480.8179 


2005, 


.3000 


+32, 


,11 


0, 


,41 


+0, 


.03 


0, 


.972 


53482.6947 


2005, 


.3051 


+33, 


,57 


0, 


,44 


-0, 


.26 


0, 


.976 


53487.7003 


2005, 


.3188 


+38, 


,96 


0, 


,44 


+0, 


.07 


0, 


.984 


53510.6947 


2005, 


.3818 


+54, 


.52 


0, 


.43 


+0, 


.13 


0, 


.025 



''Velocity uncertainties include the scale factor described in the text. 
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Table 2. Orbital solutions for HIP 50796. 



Parameter 



Spectroscopic only 



Combined" 



Constrained'^'*' 



Adjusted quantities 



P (days) 


570.70 


± 


0.73 


570.95 


± 


0.52 


570.98 


± 


0.52 


7 (kms~^) 


+24.77 




0.14 


+24.86 


± 


0.11 


+24.87 


± 


0.11 


K (kms-l) 


20.62 




0.28 


20.75 


± 


0.20 


20.76 


± 


0.20 


e 


0.6103 


± 


0.0071 


0.6110 


± 


0.0052 


0.6110 


± 


0.0051 


t^i (deg) 


313.75 


± 


0.95 


314.01 


± 


0.76 


314.02 


± 


0.76 


T (HJD-2,400,000) 


52355.2 


± 


1.4 


52355.0 


± 


1.1 


52355.0 


± 


1.1 


i (deg) 








85 


± 


13 


83 


± 


15 


f2 (deg) 








175.0 


± 


8.8 


179 


± 


10 


Aa* (mas) 








+3.7 


± 


4.0 


+3.1 


± 


4.2 


A<5 (mas) 








-8.7 


± 


1.8 


-7.7 


± 


1.7 


A/i* (mas yr-i) 








+4.6 


± 


1.4 


+4.5 


± 


1.4 


Ans (mas yr-^) 








-3.4 


± 


1.3 


-3.2 


± 


1.3 


Att (mas) 








-9.9 


± 


1.8 


-8.8 


± 


1.9 


Derived quantities 




















f(M) (Mq) 


0.258 


± 


0.011 














M2 sin j/(Mi + M2)2/3 (M0) . . . 


0.6364 


± 


0.0092 














ai sini (10® km) 


128.2 


± 


1.9 














ai (mas) 








16.9 


± 


1.5 


17.9 


± 


1.6 


(iphot (mas) 








16.9 


± 


1.5 


14.9 


± 


1.3 


a (AU)'= 








1.580 


± 


0.029 


1.582 


± 


0.032 


a (mas)*^ 








30.8 


± 


3.2 


32.5 


± 


3.3 


Ma (™as yr-i) 








-75.4 


± 


1.4 


-75.5 


± 


1.4 


//5 (mas yr-i) 








-25.0 


± 


1.3 


-24.8 


± 


1.3 


TT (mas) 








19.5 


± 


1.8 


20.6 


± 


1.9 



"■Incorporates radial velocities, HIPPARCOS measurements (abscissae residuals), and the Tycho-2 
proper motions. 

''Assumes the light from the secondary is not negligible (a brightness difference of AV = 2.5 mag); see 
text. 



'^Assumes a value for the primary mass of Mi = 0.73 Mq. 
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Table 3. HIPPARCOS abscissae residuals for HIP 50796. 



HJD 
(2,400,000+) 


Year 


V 

(mas) 


(mas) 


Cons.^ 


Corr. 
Coef.b 


O-C 
(mas) 


Orbital 
Phase 


47884.3184 


1989 


.9776 


—6. 


,60 


4, 


,60 


F 


0. 


,762 


—5. 


,24 


0, 


,170 


47884.1723 


1 OQO 

iyoy 


OTTO 


Q 

— O. 


,oo 


0. 


A n 
,4U 


N 


n 
U. 


. iDZ 


— 1 . 


Qn 


U. 


1 TC\ 

, i ; U 


47901.6313 


1990 


.0250 


+ 1. 


,61 


5. 


.Dl 


F 


0. 


.531 


—6. 


,77 


0. 


,200 


47901.5947 


1990 


.0249 


+5. 


.55 


A 

4, 


.55 


N 


0. 


.531 


—2, 


,84 


0. 


,200 


4o020.d010 


1990 


.3499 


+3. 


.81 


A 

4, 


.31 


F 


0. 


.736 


+0. 


,01 


0, 


,4Uo 


48020.1914 


iyyu 


.o4yD 


— Z. 


■ (K) 


o. 


.Oi 


N 


U. 


■ 1 oD 


— D. 


AC\ 

,4y 


n 
U, 


,4U<5 


48046.1241 


1990 


.4zUD 


+20. 


,01 


8. 


.OI 


F 


0, 


.868 


1 1 o 

+12. 


,7o 


0. 


A CQ 

,45o 


48046.0146 


1 onn 

iyyu 


A onQ 


1 on 


ni 
Ui 


1 n 
iU. 


T^ 

■ ( O 


N 


n 
U. 


.OUO 


1 1 o 

+ iz, 


, IO 


U. 


A f^Q 


48046.3068 


1990 


.4211 


+ 1, 


.86 


4, 


,72 


F 


0. 


.773 


—5. 


.37 


0. 


,454 


48046.3798 


1 onn 

iyyu 


A 01 Q 


1 Q 


tiT 


r: 
O. 


CO 

,oZ 


N 


U. 


. 1 iO 


+1, 


A A 
,44 


U. 


AKA 
,4o4 


48069.0519 


1 nnn 

iyyu 


AQQA 

.4oo4 


+0. 


.Do 


Q 
O. 


QA 

,o4 


F 


U, 


QTQ 
.O / O 


1 n 
+U. 


OK 
,20 


U. 


AC\A 

,4y4 


48069.0619 


1990 


.4834 


— 0. 


1 A 

.14 


9. 


.83 


N 


0. 


.878 


—6. 


.OO 


0. 


A(\A 

,4y4 


48069.2445 


1990, 


.4i5oy 


—2. 


on 


D. 


cn 


F 


0. 


nnn 
UUU 


—8. 


.62 


u. 


A(\A 

,4y4 


48203.1086 


1 non 

iyyu 


.?5oU4 


— i. 


on 

.c5U 


Q 

O. 


,y t 


F 


U, 


■ ( oy 


1 o 

+Z, 
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0, 
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Table 3 — Continued 



HJD 




V 






Corr. 


O-C 


Orbital 


(2,400,000+) 


Year 


(mas) 


(mas) 


Cons.'' 


Coef.t' 


(mas) 


Phase 


48760.4436 


1992.3763 


+9.38 


9.65 


F 


0.751 


+ 10.30 


0.704 


48760.4070 


1992.3762 


-2.08 


13.41 


N 


0.751 


-1.18 


0.704 


48808.2183 


1992.5071 


+10.67 


11.83 


F 


0.875 


+8.18 


0.788 


48809.0218 


1992.5093 


+4.26 


13.06 


N 


0.875 


+1.84 


0.789 


48809.2410 


1992.5099 


+3.03 


5.48 


F 


0.505 


+0.89 


0.790 


48809.0949 


1992.5095 


+1.71 


9.46 


N 


0.505 


-0.42 


0.790 


48942.7033 


1992.8753 


-2.01 


6.43 


F 


0.777 


-9.08 


0.024 


48942.7033 


1992.8753 


+2.19 


7.92 


N 


0.777 


-4.88 


0.024 



'■Data reduction consortium responsible for the measurement: F = FAST, N = NDAC (see 
ESA 1997). 

''Correlation coefficient between the FAST and NDAC abscissae taken on the same great 
circle. 
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Table 4. Parallax and proper motion components of HIP 50796. 









TT 


Source 


(mas yr""^) 


(mas yr~^) 


(mas) 


HIPPARCOS catalog (ESA 1997) . . . 


-80.0 ± 2.4 


-21.6 ± 1.9 


29.4 ± 2.7 


Initial fit (this paper) ^ 


-76.1 ± 2.2 


-24.6 ± 1.8 


19.4 ± 1.9 


iyclio-2 catalog (H0g et al. 2000) . . . 


-74.9 ± 1.8 


-25.3 ± 1.8 




Combined fit (this paper)^ 


-75.4 ± 1.4 


-25.0 ± 1.3 


19.5 ± 1.8 


Constrained fit (this paper) 


-75.5 ± 1.4 


-24.8 ± 1.3 


20.6 ± 1.9 



'^Incorporates radial velocities and HIPPARCOS abscissae residuals, but not the 
Tycho-2 proper motions. 

"^Includes the Tycho-2 proper motions, and assumes the companion contributes no 
hght. 

^Accounts for the small light contribution from the companion. 
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Table 5. Photometry for HIP 50796. 





Mv 


V -J 


V -H 


V -K 


Source 


(mag) 


(mag) 


(mag) 


(mag) 


Observed (tt from combined solution) 


7.25 ± 0.20 


2.32 ± 0.03 


2.99 ± 0.04 


3.10 ± 0.04 


Girardi et al. (2000) isochrones'* 


7.25 


1.96 


2.60 


2.67 


Henry & McCarthy (1993)=' 


7.22 


2.15 


2.70 


2.85 


Observed (tt from constrained solution) . . . 


7.37 ± 0.20 


2.32 ± 0.03 


2.99 ± 0.04 


3.10 ± 0.04 


Model for triple system 


7.14 


2.34 


2.95 


3.12 



^Predicted photometry for a single star with Mi = 0.73 Mq. 



